Mott correlated states in the underdoped two-dimensional Hubbard model: 
variational Monte Carlo versus a dynamical cluster approximation. 
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We investigate the properties of the frustrated underdoped Hubbard model on the square lattice 
using two complementary approaches, the dynamical cluster extension of dynamical mean field 
theory, and variational Monte Carlo simulations of Gutzwiller-Jastrow wavefunctions with backfiow 
corrections. We find good agreement, apart from the location of the Mott-Hubbard transition 
which differs. At small dopings, we observe a rapid crossover from a weakly correlated metal at 
low interaction strength U to a non-Fermi liquid correlated state with strong local spin correlations. 
Furthermore, we investigate the stability of the correlated state against phase separation. We 
observe phase separation only for large values of U or very large frustration. No phase separation 
is present for the parameter range relevant for the cuprates. 

PACS numbers: 71.10.Fd,71.27.+a,71.30.+h,71.10.Hf 
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I. INTRODUCTION 

The Hubbard model on the square lattice is a minimal 
model for describing electronic correlation. This model 
has played a central role for the study of high-Tc super- 
conductivity, since it is believed to capture the essential 
physics of the copper-oxygen planes in cuprate materi- 
als. Furthermore, the properties of the Hubbard model 
in the underdoped regime, i.e. in the proximity of the 
Mott insulating state at half-filling, is a stimulating re- 
search area due to the challenge in describing the physics 
of the correlated pseudogap state and its non-Fermi liq- 
uid behavior. Significant achievements in this field have 
been made, for instance, by the cluster extensions of dy- 
namical mean-field theory^—, by means of a phenomcno- 
logical theory 7 - and of the Gutzwiller approximation^, as 
well as on the basis of increasingly accurate variational 
wave functions*^. 

Recently, it has been proposed that the pseudogap and 
the superconducting phases present at finite dopingiir— 
and for onsite U values corresponding to a Mott insulat- 
ing state at half filling, can be continuously connected to 
a pseudogap and a superconducting phase at half-filling 
for U values lower than the critical U c corresponding 
to the Mott-Hubbard metal-insulator transition (MIT). 
This proposal contrasts with Anderson's concept of su- 
perconductivity as a state emerging out of a Resonant Va- 
lence Bond (RVB) stated, a prototypical Mott-insulating 
state at half filling in the absence of magnetic order. 

Another interesting feature of a correlated electron 
state in the underdoped regime is the tendency of the 
system to phase separate into an undoped state with 
strong antiferromagnetic correlations and a hole-doped 
region. Indeed, if the Hubbard model would be unstable 
to phase separation, its validity as a model to properly 
describe the development of superconductivity could be 
questioned. Phase separation occurs when the stabil- 
ity condition d 2 E(n)/dn 2 > is violated, i.e. when the 
ground-state energy E(n), as a function of electronic den- 



sity n, is not any more convex. As introduced by Emery 
et ali^, phase separation can be studied by looking at 
the hole energy Eg(6), defined as, 



Eg (6) 



E(S) - E(0) 



(1) 



where S = 1 — n is the hole density. If the hole energy has 
a minimum at a critical doping 6 C , the system is unstable 
to phase separation for S < 6 C . A difficulty when using 
the Emery criterion ([T]) is its strong dependence on the 
accuracy in the estimate of the ground state energy, as 
discussed in Ref. Il6l In particular, less accurate ground 
state energies tend to overestimate the critical 6 C below 
which the system exhibits phase separation. This is a 
consequence of the fact that good estimates for the en- 
ergy of low-doping states are hard to obtain because of 
the strong local correlations. Different approaches have 
led in the past to contradictory results for the Hubbard 
model in the parameter range relevant for hole-doped 
cuprates, ranging from absence of phase separation to 
a range of estimates for phase separation up to 10% dop- 
ingi 7 . - — . For the t — J model2i~— , i.e. for the large- U 
limit of the Hubbard model, phase separation is present 
at all dopings for large values of J, but the possible oc- 
currence of phase separation close to half filling for small 
values of J has been a long-standing debate^—. 

In this work, we make use of variational Monte Carlo 
(VMC) simulations, including backfiow correlations^, 
and the dynamical cluster approximation (DCA)22i2£ to 
address the properties of the Hubbard model on the 
square lattice with nearest and next-nearest neighbor 
hoppings in the underdoped regime both at zero and fi- 
nite temperature. As already pointed out in Ref. [ToL we 
distinguish a weakly correlated metal at small U (i.e. 
continuously connected to a half-filled metallic regime) 
and a strongly correlated state at intermediate to large 
U (i.e. when the half-filled case is insulating). We find 
significant differences between the two regimes in the 
static and in the dynamical spin correlations as well as in 
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the low-energy part of the self-energy. The above obser- 
vations obtained from two complementary approaches, 
VMC and DCA, and the violation of the Luttinger sum 
rule for the strongly correlated state reported in a pre- 
vious study 10 , hint to a non- Fermi liquid nature of the 
strongly correlated state. The two states are separated 
by a crossover line at the critical interaction U c , that 
may evolve into a first-order transition^ when doping is 
vanishing. 

Moreover, we investigate the occurrence of phase sep- 
aration in the underdoped Hubbard model. We find no 
tendency to phase separation when the Coulomb repul- 
sion U is slightly above the critical U c , regardless of 
the ratio t /t between next-nearest and nearest neigh- 
bor hoppings, t' and t respectively. Our results there- 
fore indicate that for intermediate values of the Coulomb 
repulsion (U/t ~ 6 — 8), the Hubbard model is stable 
against phase separation, irrespective of the value of t' jt. 
The actual value of t'/t is sensitive to the details of the 
high-T c cuprate compound investigated and can be calcu- 
lated, for example, in density functional theory, leading 
to t'/t ~ [—0.1, —0.4]. A larger degree of frustration has 
been shown to correlate, in general, with higher critical 
temperatures^. 

If we increase the electron repulsion U at constant t'/t, 
the system is, in contrast, found to be unstable towards 
phase separation. However, since increasing the varia- 
tional accuracy tends to suppress phase separation, we 
cannot exclude that an even more accurate approach 
could eliminate or significantly reduce the tendency to 
phase separate also for larger values of the Coulomb re- 
pulsion. 

The paper is organized as follows: in Section [Hi we 
introduce the Hamiltonian and describe the VMC and 
the DCA approaches; in Section Hill we compare ground- 
state energy and the density of double occupancies within 
the two methods, in Section ITVl we discuss the non-Fermi 
liquid properties of the underdoped region at U > U c , in 
contrast with the weakly correlated metal at U < U c . In 
Section [V] the occurrence of phase separation is discussed 
and finally we present our conclusions in Section IVI1 

II. MODEL AND METHODS 



We consider the frustrated Hubbard model with ex- 
tended hopping on a two-dimensional square lattice, 

U = -tjy ia c jc - t' J2 <L c j„ + Hc - + U^2n i tn ii , 
(ij)o 1 

(2) 

where c| (c ia ) denotes the electron creation (annihila- 
tion) operator of one electron on site i with spin a =~f, J,, 
(ij) and ((ij)) indicate nearest and next-nearest neighbor 
sites respectively; n ia = c\ a c ia is the electron density; t 
and t' are the nearest and next-nearest neighbor hopping 
amplitudes, and U is the on-site Coulomb repulsion. 
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FIG. 1. (Color online) (a) Double occupancy D and (b) en- 
ergy as a function of electronic density n for U/t = 7.0 and 
t' — — OAt. Results are obtained using a VMC approach on an 
L = 162 lattice size (blue circles) and within DCA by means 
of a 2 x 2 plaquette in fc-space for a temperature T/t = 0.05 
(red squares). Errors in VMC are smaller than the symbol 
size. Lines are just guides to the eye. 

A. Variational Monte Carlo 

Variational Monte Carlo (VMC) simulations allow to 
perform non-perturbative calculations at zero tempera- 
ture for one- and two-dimensional correlated and frus- 
trated systems. VMC simulations are based on numer- 
ically sampling expectation values over a variational es- 
timate of the ground-state wave function. Here we will 
use a powerful variational ansatz for frustrated electron 
systems which has been strictly tested by comparing 
extensively to analytical and numerical exact limiting 
cases. Our variational ansatz has been tested in particu- 
lar against Bethe ansatz predictions for the Luttinger liq- 
uid exponents in one dimension^, with respect to Lanc- 
zos results, in terms of wave function overlap, for a two- 
dimensional 18 site cluster—, and with respect to the 
results of density matrix renormalization group studies, 
in terms of the ground-state energy, for the ID J\ — J2 
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FIG. 2. (Color online) (a) Double occupancy D and (b) en- 
ergy as a function of U/t for electron density n = 0.95 and 
t' — — OAt. Data are obtained by a VMC approach on an 
L — 162 lattice size (blue circles) and within DCA by means 
of a 2 x 2 plaquette in fc-space for a temperature T/t = 0.05 
(red squares). Errors in VMC are smaller than the symbol 
size. Lines are just guides to the eye. 



mode! 28 . 

The variational ansatz consists of three components. 
In a first step, we construct uncorrelated non-magnetic 
wave functions given by the ground state |BCS) of a su- 
perconducting Bardeen-Cooper-Schrieffer (BCS) Hamil- 
toniai>2i-J£: 



Hbcs = & c L c k<r + AkC l 



(3) 



where both the free-band dispersion and the pairing 
amplitudes Afe are variational functions. We use the 
parametrization 

£fc = — 2i(cos k x + cos k y ) — At' cos k x cos k y — \l (4) 
A fe = 2A B cs(cosfc K - cosk y ), (5) 

where the effective hopping amplitude r, the effective 
chemical potential fi, and the local pairing field Abcs 



are variational parameters to be optimized in order to 
minimize the variational energy. The parameter t is kept 
fixed to set the energy scale. We point out that large 
electronic correlations, as in proximity of a Mott insula- 
tor, lead to a strong renormalization of t' with respect to 
the bare Hamiltonian value t' , see Ref . |T(J 

The correlated state I'I'bcs), without backflow terms, 
is then given by |*bcs) = J|BCS), where J = 
exp(— 1/2 y~] i: j VijUiUj) is a density-density Jastrow factor 
(including the on-site Gutzwiller term vu), with the uy's 
being optimized for every independent distance \i — j\. 
Notably, within this kind of wave function, it is possi- 
ble to obtain a pure (i.e., non-magnetic) Mott insula- 
tor for a sufficiently singular Jastrow factor v q ~ l/q 2 
(where v q is the Fourier transform of ) , while a super- 
conducting (metallic) state is found whenever v q ~ l/q 
and Abcs > (A BC s = 0)24. 

A size-consistent and efficient way to further improve 
the correlated state |^bcs) I0r large on-site interactions 
is based on backflow correlations. In this approach, each 
orbital that defines the unprojected state |BCS) is taken 
to depend upon the many-body configuration, in order 
to incorporate virtual hopping processes^, in particular 
the recombination of neighboring charge fluctuations is 
favored. This is a substantial improvement with respect 
to Jastrow factors, where electron-electron correlation is 
included only via a multiplicative term. All results pre- 
sented here are obtained by fully incorporating the back- 
flow corrections and optimizing individually^ every vari- 
ational parameter in and Afe, in the Jastrow factor J , 
as well as for the backflow corrections. Calculations are 
performed on 45° tilted clusters with L — 162 lattice sites 
and periodic boundary conditions. 



B. Dynamical Cluster Approximation 

DCA is the cluster extension of single-site dynamical 
mean field theory (DMFT) 38 , which includes, to a certain 
degree, momentum dependencies. Since the hopping ma- 
trix for sites within the considered cluster and for sites on 
different clusters is the same, in contrast to the cellular 
DMFT approach^, the DCA self-consistent equation can 
be written in momentum space with the assumption that 
the self-energy is constant in the Brillouin zone sectors 
that are considered. The DCA self-consistency equation 
is given as 



G CT (K, iu n ) 



N ^ iu 

K 



1 



M~ e K+K _ £<t(K,MJ„)' 

where N is the number of k points in each Brillouin zone 
sector, [i the chemical potential, K is the cluster momen- 
tum, e K , g the dispersion relation, u n are the fermionic 

Matsubara frequencies, and where the summation over K 
is performed in each Brillouin zone sector. In our calcula- 
tions, we employed the DCA cluster with N c = 4, where 
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FIG. 3. (Color online) Evolution of the static structure fac- 
tor correlations S(q) = (s- q s q ) at Q = (n, n) as a function 
of U/t for n = 0.975 (blue squares) and n = 0.95 (red cir- 
cles) at t'/t = —0.4. The appearance of the non-FL region is 
characterized by a rapid increase in the short-range antiferro- 
magnetic correlations. Data are obtained by means of VMC 
simulations on an L = 162 lattice size. 

K = (0,0), (0,7r), (tt, 0), and (w,tt). The converged self- 
energy S(t(K, iuj n ) is evaluated by means of Eq. © and 
the Dyson equation. We employ the interaction expan- 
sion continuous-time quantum Monte Carlo approach as 
an impurity solve r 40 ' 41 . All calculations presented here 
are for a temperature T/t = 0.05 and we perform more 
than 10 7 QMC samplings for the impurity Green's func- 
tion G(iuj n ), in order to keep the QMC statistical errors 
smaller than 5 x 10~ 3 for the first Matsubara frequency. 



while data at strong coupling (U/t ~ 7.0) and at weak 
coupling (U/t < 3.5) are in good agreement. This dis- 
crepancy is related to the fact that VMC and DCA pre- 
dict different values for the location of U c /t of the metal- 
insulator transition at half-filling. The singlet RVB state 
is dominant in the DCA on a 2 x 2 cluster, thus favor- 
ing correlated states, and also the VMC results can be 
slightly dependent on the accuracy of the trial wave func- 
tion. The metal-insulator transition at half-filling, which 
is weakly first-order, takes place at U c /t — 5.8±0.2 within 
VMG±°., while for DCA, on a 2 x 2 cluster, the critical 
U c /t is approximately estimated in U c /t = 4.6 ± 0. 4 12 ' 42 . 
Note that the metal-insulator transition in VMC has 
been determined by looking at the static structure fac- 
tor N(q) = (n- q n q ), where n q = 1/ 'VLJ2 r ,a eiqrn r,<? is 
the Fourier transform of the particle density. Indeed, 
the metallic phase is characterized by N(q) ~ q for 
q — > 0, which implies a vanishing gap for particle-hole 
excitations. On the contrary, in the insulating phase, 
N(q) ~ q 2 for q — > 0, implying the fact that the charge 
gap is finite. 

The key point of Fig. [2] is however the fact that signa- 
tures of the metal-insulator transition occurring at half 
filling are strongly visible at low doping for both ap- 
proaches, independently of the actual location of the MIT 
predicted in each method. The curves for both the double 
occupancy D and the energy show kinks for interactions 
U close to the respective critical U c . We thus infer that 
a rapid crossover between a weak-coupling Fermi-liquid 
(FL) and an intermediate to large U non-FL regime oc- 
curs at small but finite dopings. This conclusion is consis- 
tent with the observed behavior of the one-particle self- 
energy and of the two-particle correlation functions which 
we will discuss next. 



III. ENERGY AND DOUBLE OCCUPANCIES 



In Fig. [TJ we compare the average number of double 
occupancies per site D = (n^fn^j.) and the energy, for 
U/t = 7.0 and t' = — OAt as a function of the electron 
density n, as obtained from VMC and DCA simulations. 
For the DCA calculations, the energy is calculated by 

E = I? E [efcG ff (fc,^„)] e ^" 0+ + UD, (7) 

where the index k runs over the first Brillouin zone 
and we considered the asymptotic behavior of the self- 
energy in the limit of large Matsubara frequencies iuj n : 
S<j(K, iuj n ) U 2 n a (l — n a )/i(jj n . Both data sets are 
in good agreement, suggesting that both approaches al- 
low for a reasonably accurate description of the two- 
dimensional model. 

In Fig. [2l we present the double occupancy D and the 
energy close to half filling, for n — 0.95 and t' = —OAt, 
as a function of the interaction strength U/t. We ob- 
serve that the results differ in the range U/t = 4 — 6, 



IV. NON-FERMI LIQUID VS. FERMI LIQUID 
PROPERTIES 



In this Section, we explore in more detail the FL and 
non-FL properties of the frustrated Hubbard model close 
to half filling, using both the VMC and DCA approaches. 
Long-range static correlations can be evaluated within 
VMC but not within DCA, due to cluster-size restric- 
tions. Dynamical quantities like the one-particle self en- 
ergy can however be calculated within DCA and not with 
VMC; the two approaches complement each other nicely. 

The nature of the non-Fermi liquid region has been 
characterized, using VMC simulations^ ., by a strong 
renormalization of the underlying Fermi surface and a 
small violation of the Luttinger sum rule. Here, we as- 
sess the magnetic properties at low dopings via the static 
structure factor, defined as 

S(q) = jJ2 eiq(Rm ~ Rn) ( S ™ S -)> ( 8 ) 

m,n 

where is the z-component of the spin operator on site 
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FIG. 4. (Color online) Upper panel: The imaginary part of 
the lowest (red circles) and the second lowest (blue squares) 
self-energy values in Matsubara frequencies as a function of 
U/t at T/t = 0.05 with n = 0.95 and t' = -OAt. Data have 
been obtained by DCA using a 2 x 2 cluster. The crossing 
point of both self-energies occurs around U/t = 4.5, which is 
the same critical U c /t where the kink in the double occupancy 
data of Fig.[2]is located. The quantity a is defined in Eq. 
Lower panel: The imaginary part of the self energy as a func- 
tion of the Matsubara frequencies at U/t — 3.6 and n = 0.94 
(FL, red squares) and at U/t = 5.0 and n = 0.96 (non-FL, 
blue circles). Data have been obtained by DCA using a 2 x 2 
cluster. 



m and where L denotes the total number of sites. The 
presence of (short-range) antiferromagnetic correlations 
is signaled by the appearance of a (non-diverging) peak 
in S(q), located at Q = (ir, ir). 

As shown in Fig. [31 for the two electron densities 
n = 0.975 and n = 0.95 at t'/t = -0.4, the non-Fermi 
liquid state at U/t > 6 is characterized by antiferro- 
magnetic correlations which are substantially enhanced 
with respect to the weakly correlated metallic phase at 
U/t < 6. For the smaller doping, n — 0.975, the two 
regimes are clearly separated by a rapid increase in the 
value of S(Q), which could be compatible with a first- 
order transition, while for the larger doping the jump is 
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FIG. 5. (Color online) The local dynamical spin suscepti- 
bility Im(x sp H) for U/t = 3.6 and n = 0.94 (FL, open 
red squares), U/t = 5.0 and n — 0.96 (non-FL, filled green 
squares) and U/t — 7.0 and n = f.O (Mott insulator, open 
blue circles), as a function of the real frequency u/t at 
T/t = 0.05 and t! = —OAt. A Pade approximation has been 
employed for the analytical continuation of X sp ( T )- 



less evident and the observed rise in the strength of the 
short-ranged spin-spin correlations is more in agreement 
with a smooth crossover. 

Next, we plot the imaginary part of the lowest and the 
second lowest Matsubara frequency self-energy values as 
a function of U/t at T/t = 0.05, n = 0.95 and t' = -OAt 
(see Fig.Sl upper panel), obtained using DCA. We define 



Im(S(ia>i)) — Im(S(icJo)) 



wi - UJQ 



(9) 



where Im(£(m>o)) is the imaginary part of the lowest 
Matsubara frequency self-energy value and Im(E(io;i)) is 
the imaginary part of the second lowest Matsubara self- 
energy value. Negative ratios a < indicate (quasi)-FL 
behavior, while positive ratios a > suggest a non-Fermi 
liquid state. Indeed, the value of the imaginary part of 
the self-energy converges to zero (or to small values due 
to the effect of temperature) for iuj n — > in the Fermi- 
liquid regime, while it is monotonically decreasing for 
iuj n — > in the non- Fermi liquid state, see Fig. 01 lower 
panel. From the data presented in Fig. 21 we find a criti- 
cal U c /t = 4.5 between FL and non-FL states, which is in 
agreement with the critical U c extracted from the small 
kink in the double occupancy D, compare Fig. [2] 

Finally, we present in Fig. [S] the local dynamical spin 
susceptibility Im(x sp (u;)) obtained within DCA with 10 8 
QMC samplings, by performing the analytical continua- 
tion of 



X sp (t) = (S z (t)S z (0)), 



(10) 



where S z = \{n^ — njj and r is the imaginary time, with 
the help of the Pade approximation. The Im(x sp (w)) 
around u>/t — 0.5 is strongly suppressed in the FL region 



6 



-1.4 



q- -1.8 
UT 

{$L -2 
LU 

-2.2 
-2.4 



III 


- 
- 


^^^^ 

♦ — -^^r'^ 


t 


: 


(y/r=5« : 




U/f=6» - 




t//f=7T . 


, i . . i . . i . 


U/t=10* - 

III" 



0.03 0.06 0.09 0.12 0.15 

5 



Uj 



-1.6 



-1.8 



-2.2 



-2.6 



1 1 1 1 1 

- 


i 1 1 i 1 1 i 1 1 


- 






^^r/t=ow - 




t'/t=-0.2 ♦ - 




t'/t=-0.4» - 


■ i i i i 


t'/t=-0.75w 1 
i , , i , , i , , 



0.03 



0.06 0.09 

5 



0.12 0.15 



FIG. 6. (Color online) The hole energy E S (S) = (E(6) - 
E(0))/8 as a function of doping for different values of U/t at 
t'/t = —0.4. Data are obtained by means of VMC simulations 
on an L — 162 lattice size 



at U/t = 3.6 and n = 0.94. This low-frequency peak is 
dominant in the non-FL region at U/t = 5.0 and n = 0.96 
and in the Mott insulating region at U/t — 7.0 and n = 1. 
We relate the enhancement of the low-energy peak in 
non-FL and Mott insulator regions to the formation of 
short-range antiferromagnetic correlations, as presented 
in Fig. [3J in terms of increased RVB-type singlet pair- 
ing. These results demonstrate the complementarity of 
the DCA and VMC methods in identifying the region of 
possible non-FL behavior and they are compatible with 
studies of the Kagome and the triangular lattic o 43 ' 44 . 



PHASE SEPARATION 



We investigate now the possible occurrence of phase 
separation in the hole-doped regime n < 1, by consid- 
ering the hole energy defined in Eq. ((T|). The system is 
unstable to phase separation for 5 < S c , if the hole energy 
has a minimum at a critical doping S c , while a monoton- 
ically increasing hole energy corresponds, on the other 
hand, to an energetically favorable homogeneous solu- 
tion for the doping levels investigated. Physically, phase 
separation is driven by magnetic correlations, which can 
be substantially increased, in the phase without holes, at 
the expense of the kinetic energy. 

In Fig. [HI we focus on the case t'/t — —0.4, and vary 
the interaction strength U/t. As expected, we find that 
no phase separation is possible when the half-filled case is 
metallic (U/t — 5), due to the absence of a relevant mag- 
netic energy scale. Phase separation is also not observed, 
interestingly, when U is only slightly above the critical 
interaction U c (U/t = 6, in Fig. [6]), at least for small 
but finite doping levels 6 > 0.02. At very low dopings 
S < 0.02 (and U/t — 7) phase separation could possi- 
bly be present. However, since the model is expected to 



FIG. 7. (Color online) The hole energy E S (S) = (E(S) - 
E(0))/S as a function of doping for the four cases U/t = 
6, t'/t = (blue squares), U/t = 6,t'/t = -0.2 (pink dia- 
monds), U/t = 6, t'/t = -0.4 (red circles) and U/t = S,t'/t = 
—0.75 (black triangles). The data has been obtained by means 
of VMC simulations on an L — 162 lattice. 



exhibit a magnetic instability close to half-filling, a pos- 
sible phase separation occurring for 6 < 0.02 would be 
masked by long-range magnetic order. For larger values 
of U/t phase separation seems however to be energeti- 
cally favorable for a wider range S < 0.08 of doping lev- 
els. The data for larger electron repulsion (U/t ~ 20) are 
not shown in Fig. O but they exhibit trends similar to 
the case U/t = 10. However, we cannot exclude that bet- 
ter estimates of the ground-state energy can reduce this 
tendency to phase separate, as detailed in Ref. Il6l This 
is a consequence of the fact that good estimates for the 
energy of low-doping states are hard to obtain because 
of the strong local correlations. 

In Fig. \7\ we show the hole energy as a function of 
doping for four values of the next-nearest neighbor hop- 
ping, ranging from t'/t — to t'/t = —0.75. The value 
of the Coulomb repulsion is slightly above the critical 
U c , that is located between U/t — 5 and U/t — 6 for 
-0.4 < t'/t < 0, while it is located between U/t = 7 
and U/t = 8 at t'/t = —0.75. In all the cases relevant 
for the cuprates, t'/t ~ [—0.1, —0.4], no phase separation 
occurs. Only in a small range 5 < 0.02 a tendency for 
phase separation is observed for t'/t — —0.75. Thus, our 
data suggest that for the cuprates the Hubbard model 
is not unstable against phase separation when the value 
of the electronic repulsion is chosen to be close to the 
critical interaction U r of the Mott-Hubbard transition. 



VI. CONCLUSIONS 



When investigating correlated electron systems numer- 
ical or analytical approximations are generically neces- 
sary and the accuracy of the employed approach is no- 
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toriously difficult to control. Here we compare results 
obtained by two complementary approaches, DCA and 
VMC. We find good overall agreement, apart from the 
value of the critical Hubbard- U for the Mott-Hubbard 
metal-insulator transition at half filling. We use the com- 
plementary information, static long-range correlations 
provided by VMC and dynamical properties provided by 
DCA, to investigate the physics at finite but low doping 
levels. We find that the crossover from a weakly corre- 
lated electron state at intermediate to small values of U 
to a non-Fermi liquid state at intermediate to large val- 
ues of U is characterized by a strong increase in local 
magnetic correlations. In this respect, we do not find 
evidence for non-Fermi liquid properties below the crit- 
ical Hubbard-JJ, as suggested instead by a recent DCA 
stud y 11 ' 12 . Our result is a solid feature of the VMC ap- 
proach— The application of DCA in the present work 
has been oriented to supplement the dynamical proper- 
ties that are missing in VMC and we do not exclude 
that other studies within DCA can show evidence of non- 



Fermi liquid properties also below the critical U, though 
this is in contradiction with VMC. This remains a con- 
troversial issu o 12 ' 45 and needs further study. In Sec. [Vj 
we investigate the stability of the non-Fermi liquid state 
against phase separation and find it to be stable for all pa- 
rameters relevant for the cuprates. These studies further 
prove the synergies obtainable when using complemen- 
tary methods for the study of frustrated and correlated 
electron systems. 
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